Origin of thermoelectric response fluctuations in single-molecule junctions 



Yonatan Dubi 1 '^ 

1 Department of Chemistry and the Use Katz center for Nano-Science, 
Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel 
(Dated: November 21, 2012) 

The thermoelectric response of molecular junctions exhibits large fluctuations, as observed in 
recent experiments [e.g. Malen J. A. et ai., Nano Lett. 10, 3406 (2009)]. These were attributed 
to fluctuations in the energy alignment between the highest occupied molecular orbital (HOMO) 
and the Fermi level at the electrodes. By analyzing these fluctuations assuming resonant transport 
through the HOMO level, we demonstrate that fluctuations in the HOMO level alone cannot account 
for the observed fluctuations in the thermopower, and that the thermo- voltage distributions obtained 
using the most common method, the Non-equilibrium Green's function method, are qualitatively 
different than those observed experimentally. We argue that this inconsistency between the theory 
and experiment is due to the level broadening, which is inherently built-in to the method, and 
smears out any variations of the transmission on energy scales smaller than the level broadening. 
We show that although this smmearing only weakly affects the transmission, it has a large effect 
on the calculated thermopower. Using the theory of open quantum systems we account for both 
the magnitude of the variations and the qualitative form of the distributions, and show that they 
arise not only from variations in the HOMO-Fermi level offset, but also from variations of the local 
density of states at the contact point between the molecule and the electrode. 



I. INTRODUCTION 

Improving the thermoelectric (TE) energy conversion 
efficiency relies on increasing the thermoelectric response 
(Seebeck coefficient) and reducing the thermal conduc- 
tance of the device^. Metal-single molecule-metal junc- 
tions ('molecular junctions') seem promising in both as- 
pects: their thermal conductance is small due to the 
mismatch between the vibrations of the electrodes and 
the molecule, and their TE response should be large 
due to the well-defined resonant structure of the elec- 
tron transport (via the molecular HOMO or LUMO 
levels) 3 . This observation, along with the notion that 
TE response can shed light on transport mechanisms in 
molecular junctioniP^, have initiated in large number of 
experimental^ anc i theoretical studies^HSl on TE ef- 
fects in molecular junctions. 

Most notable are a series of impressive experiments 
in which TE conv ersi on using a single-molecule junction 
was demonstrate d 3 ^ " 10 [ In these experiments, a single 
molecule (usually Benzene rings with various end groups) 
is trapped between an Au substrate and an Au scanning 
tunneling microscope (STM) tip, which arc held at a con- 
stant temperature difference AT. A voltage bias AU 
(here called 'thermo- voltage') is then applied between the 
tip and substrate to reach a state of zero current flowing 
through the molecule. The Seebeck coefficient is defined 
as (minus) the slope of AV(AT) as a function of AT in 
the linear response regime (i.e. AT — » 0). In the ex- 
periments the value of AU is strongly fluctuating, and 
repeating the experiment many times results in a broad 
distribution of AU, as can be seen in Fig. 3 of RefP (and 
inset of Fig. 2jb) here). Using the typical AU (where 
the voltage distribution displays a maximum), the au- 
thors of RefP^ find typical Seebeck coefficients of the 
order of ~ 10/iV/K for single molecule junctions. 



The observed variations in the thermoelectric response 
do not seem to be simple (Gaussian) noise, and the 
distributions of the thermo- voltage AU may have a 
well-defined double-peak structure^. The authors of 
Ref. [HI concluded that the fluctuations (also reported 
in thermoelectric measurements of metal-fullerene-metal 
junctions^!) are caused by variations of the position of 
the HOMO level, which are translated into variations in 
the energy offset between the HOMO level and the Fermi 
energy of the contacts, E — Ehomo~Ef, as the junction 
reconstructs. Assuming that the transport is dominated 
by a single resonant level and using the Landauer for- 
malism (or Green's function method) for thermoelectric 
transport (see, e.g. Refs. 121281 ) they quantified the vari- 
ations 8E in E and found them to be 5E ~ 2.5eU, of the 
order of the average offset E itself. 

In this paper we argue that besides the variations in 
the offset energy between the HOMO level and the Fermi 
energy, there is an additional contribution to the fluctua- 
tions in the thermopower, which is the strong variations 
of the LDOS (as a function of energy) at the point of con- 
tact between the molecule and the electrodes. We use a 
toy model to verify this effect. The main results are (I) 
Within the Landauer or Green's function formalism for 
thermoelectric transport, any variations in the transmis- 
sion function which are smaller than the self-energy (or 
level broadening) are suppressed. The result is a smooth 
transmission function and a smooth thermopower (as a 
function of energy). (II) The thermopower is extremely 
sensitive to small fluctuations of in the transmission func- 
tion. As a result, the suppression of fluctuations in the 
transmission function, inherent to the Green's function 
method, strongly affects on the calculated thermopower, 
and results in an over-estimation for the fluctuations in 
the HOMO level-Fermi level energy offset fluctuations. 
(Ill) The thcrmo-voltage histogram obtained using the 
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Green's function method is qualitatively different from 
the histograms obtained experimentally, and (IV) all of 
the above draw-backs of the Green's function formal- 
ism can be overcome by using the open quantum system 
approach to calculating the thermo-electric response in 
molecular junctions. 

The paper is organized as follows. In Sec. [IT] we derive 
the estimation for the fluctuations in the HOMO-Fermi 
level energy offset from the experimental values of ther- 
mopower for typical molecular junctions, using a simple 
description of the molecular junction as a resonant level 
and using the Landauer formula. In Sec. |III| we present 
the open quantum systems approach, and using a sim- 
ple model for the molecular junction we show that this 
approach for calculation of thermopower yields results 
which are in qualitative agreement with the experimen- 
tal results. In Sec. llVlwe show that, for the same model. 



the Green's function approach gives very different results, 
due to the inherent suppression of fluctuations built in to 
the method (due to the presence of a self-energy of level 
broadening). We demonstrate that the thermopower is 
very sensitive to fluctuations, and thus the suppression of 
(small) fluctuations in the transmission strongly affects 
the calculation of thermo-power. Sec. [V] is devoted to 
summary and discussion. 



II. ANALYSIS OF THERMOPOWER 
FLUCTUATIONS USING THE LANDAUER 
FORMULA FOR A RESONANT LEVEL 

As a starting point, we consider the data presented 
in Fig. 3 of Ref. O (also inset of Fig. [2] here), where 
the thermo-voltage histograms of a 2',5'-dimethyl-4,4"- 
tribenzenedithiol (DMTBDT) molecular junction are 
presented. We observe several notable features: (i) Noise 
is present in the measurement: even at AT w there is 
a finite thermo-voltage distribution, with an average sig- 
nal offset of a few tens of microvolts, (ii) the typical form 
of the thermo-voltage distribution seems Gaussian, (iii) 
a double-peak structure of the distribution emerges, in- 
dicating two dominant values of the Seebeck coefficient, 
Si ~ 4.6^V/K and S 2 ~ 15/iV/K. 

In order to estimate the variation in the molecular 
HOMO required to obtain these values, we use the stan- 
dard Landauer formalism . Within this formalism, the 
conductance is given by G = 2 j^T(Ep) = Gqt(Ep), and 
the Seebeck coefficient can be approximated as 
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where e is the electron energy, r(e) is the electron trans- 
mission function, e is the proton charge, kg the Boltz- 
mann constant and T is the average temperature. 

We proceed by assuming that the transport through 
molecule is characterized by a transmission through a res- 
onant level. To justify this, as well as the toy model con- 
sidered in the next section, we briefly review the litera- 



ture on theoretical calculations of thermopower in molec- 
ular junctions. The most common theoretical tool to cal- 
culate the transmission (and from it the thermopower) 
is the combination of non-equilibrium Green's function 
(NEGF) and density-functional theory (DFT)P1 where 
the transmission function is evaluated using the Green's 
function as obtained from DFT. Within this approach, 
the transmission function is calculated as a function 
of gate voltage, and from it the thermopower is calcu- 
lated, also as a function of gate voltage (see, e.g. Refs. 
II 11161 24 25 29 ) . While in transport experiments gating 
of a molecule can be achievecP^U, such a dependence of 
the thermopower on the gate voltage was never studied 
experimentally, and experiments always carry a statisti- 
cal nature as discussed above. However, to our knowledge 
the statistical nature of thermopower in molecular junc- 
tions has never been addressed (although it has been ad- 
dressed in the context of transmission^ or thermopower 
of atomic wireiP^). Nevertheless, to justify our model we 
observe that in the theoretical calculations, the transmis- 
sion function seems always to resemble a re sonant level. 
at least close to the HOMO or LUMO levels F'^ I ^ EHSl 
Thus, considering a resonant level model is a common 
description even of realistic molecules with electron in- 
teractions. 

For a given resonant molecular level with energy 
Ejjqmo an d contact Fermi energy Ep, the typical form 
of the transmission function is a resonant Lorentzian, 
T {E) = E F +r 2 , where E = E HO mo - E F is the en- 



ergy offset and T is the contact-induced effective level 
broadening. The resulting Seebeck coefficient is of the 
form 
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Noting the double-peak structure of the voltage distri- 
bution, we thus consider that there is a typical shift in 

the HOMO level 5e. Taking the conductance value to be 

r 2 

G ~ O.OlGo, one needs to solve simultaneously g2 i p2 — 
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3e E 2 +r 2 ""r* ' ' " v ' 3e 

15/xV/K . These equations yield (at ambient tempera- 
ture) T - 0.31 eV, E ~ 3.1eV and Se ~ 2.2 eV, such 
that Se/E ss 0.7, similar to the value obtained in Ref. |9] 
This means that the molecular level alignment with the 
contact Fermi level is changing by ~ 2 electron-volts at 
every reconstruction of the junction. This large variation 
(of the order of E itself) was attributed to variations in 
the junction contact geometry and intermolecular inter- 
actions. However, a recent study of transition voltage 
spectroscopy in molecular junctions (see Ref. [35] for de- 
scription of the method) demonstrated that the varia- 
tion of E are of the order of ~ leVpS much smaller than 
the variation required to generate the value above. This 
discrepancy suggests that there is an additional factor 
contributing to the variations in thermopower. 

The discrepancy describes above hints that there is a 
flaw in the analysis of the thermopower fluctuations using 
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the Landauer formula with the resonant level model. As 
we will show in the following sections, the flaw is that the 
Landauer formula inherently suppresses any variations 
in the transmission function (as a function of molecular 
level energy E) which are on an energy scale smaller than 
the level broadening (or electron self-energy), and since 
the thermopower is very sensitive to variations in the 
transmission, this suppression strongly affects the ther- 
mopower. Thus, fluctuations in E sample a thermopower 
function S(E) which is "too smooth", resulting in the 
over-estimate of SE. 



III. OPEN QUANTUM SYSTEM APPROACH 
TO THERMOPOWER 

In order to identify the origin of variations in the 
thermopower, we study a toy model of a molecular 
junction. We use the method of open quantum sys- 
tems (OQS ) (wh ich was described in detail in previous 
publication d 1 5 * 37 j) to show that the thermopower S(E) 
exhibits variations as a function of E which are on the 
same energy scale as the variations in the LDOS at 
the molecule-lead point of contact, and show that the 
thermo-voltage histogram resulting in fluctuations in E 
is qualitatively similar to the experimental ones. In the 
next section we will compare the result ontained from the 
OQS method to those obtained by using NEGF method. 

Our toy model consists of two finite electrodes with a 
molecular bridge between them (upper panel of Fig. [TJ , 
attached to reservoirs with different temperatures Tl and 
Tr for the left and right edges respectively. The molecule 
is described by a simple chain with four atomic orbitals. 
The OQS method enables one to calculate the electronic 
density in this non-equilibrium situation (a finite temper- 
ature difference between the electrodes). From the elec- 
tronic density, one can calculate (via the Poisson equa- 
tion) the voltage difference AV^ between the electrodes as 
a function of the temperature difference AT = Tr — Tl ■ 
Repeating the calculation for different values of AT we 
then find a curve Ay (AT), and the linear slope is the 
thermopower (or Seebeck coefficient) S = — . 

Essentially, OQS theory is a mapping of the system 
Hamiltonian H (which in general includes all electronic 
degrees of freedom, electron interactions, junction geom- 
etry etc.) onto a master equation for the single-particle 
density matrix of the Lindblad fornPZl, 

p = -i[H ) p]+L[p] . (3) 

We consider a non-interacting tight-binding Hamiltonian 
of the form H = % L + Hr +'Hd + T~L c , where 

H L , R =-t (4 c j+h-c) (4) 

(i,j)eL,R 

are the tight-binding Hamiltonians of the left and right 
leads respectively (t is the hopping integral, which is 



taken as the unit energy, t = leV from here on), and 

n d = -t Y ( c h + h.c) + Y( E - Ma (5) 

(i,j)ed ied 

is the Hamiltonian for the wire, which includes the usual 
hopping integral and an energy E which can be tuned 
(for instance using a gate electrode experimentally) . The 
energy is measured with respect to the electrodes' Fermi 
energy p which is set as the zero energy (/x = 0). 

The coupling between the electrodes and the molecular 
wire is described by 

"He = {9LC L C d fi + 9RC R C dl L d + h.c.) (6) 

describes the coupling between the left (right) lead to 
the wire, with c L r R \ being the creation operator for an 
electron at the point of contact between the left (right) 
lead and the wire, and Cd,o (cd,L d ) destroys an electron 
at the left-most (right-most) sites of the wire (we take 
here gL = qr = g). The external environment (s) are 
accounted for in Eq. [3]by the second term, which has the 
form 

C\p] = -\{V*Vp + pV*V) + VpV* , (7) 

where V are the Lindblad ^-operators which encode the 
properties of the environment, i.e. its temperatures and 
position (i.e. left or right electrode). An appropriate 
form for the ^-operators is^ 

V&'* = j7&*&'«>M\kM\ , (8) 

where fjj' R \ e k) = 1/ ( ex P ( RbTlr ) ~^ are ^ ae F erim 
distributions of the left and right leads (with the cor- 
responding temperature), p the chemical potential, and 

7jyfeJ are the overlap integrals between the fc and k' 
states on the left (L) and right (R) edges of the elec- 
trodes. This form for the ^-operators guarantees that at 
equilibrium (i.e. Tl = Tr) the diagonal elements of the 
density matrix are described by a Fermi function. 

We note that the Lindblad equation assumes a Markov 
approximation for the environment, which is reasonable 
since we are considering a system at room tempera- 
ture (where quantum memory effect of the environment 
should not be important), and since we are interested in 
the steady state (and not the dynamics). The numer- 
ical calculation is performed in the following way: (i) 
From the tight-binding Hamiltonian the single-particle 
states and energies are calculated, and the ^-operators 
are constructed according to Eq. Q for different elec- 
trode temperatures Tl and Tr. (ii) The ^-operators are 
inserted to Eq. (pj) and Eq. Q, which are then solved in 
the steady-state (i.e. for ,0 = 0). (iii) From the diagonal 
elements steady-state solution for p and the wave func- 
tions, the electron density is calculated, (iv) By solving 
the Poisson equation, the voltage at the center of the elec- 
trodes is calculated. The voltage difference between the 
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Figure 1: Upper panel: schematic representation of the molec- 
ular junction, (a) Seebeck coefficient S (solid blue line) and 
electron occupation of the molecular chain (dashed orange 
line) as a function energy E (see text for numerical param- 
eters). As the molecular levels cross the electrode Fermi 
energy (which is set as the zero energy) the occupancy of 
the dot change, and correspondingly the Seebeck coefficient 
changes sign. Oscillations of the Seebeck coefficients between 
the changes of sign are observed, (b) Same as (a) for dif- 
ferent values of the electrode-molecule coupling, 200, 100 and 
10 meV. The coupling has little effect on the Seebeck coeffi- 
cient, (c) temperature dependence of the Seebeck coefficient 
(at E — 0) for various molecular chain lengths, exhibiting a 
non-monotonic temperature dependence. 



electrodes is the thermo-voltage, as it is induced by the 
temperature difference, the slope of the thermo-voltage 
as a function of the temperature difference is the ther- 
mopower. 

In Fig. [l we plot the the Seebeck coefficient (solid blue 
line) and the electron occupation of the molecular chain 
(dashed orange line) of a model molecular junction as 
a function of a gate potential applied to the junction, 
i.e. as a function of the molecular energy level with re- 
spect to the Fermi level (which is set as the zero en- 
ergy) . The junction is composed of a series of four atomic 
orbitals connected to square two dimensional electrodes 
(upper panel). Numerical parameters are: electrode size 
L 2 = 25 x 25, tight-binding hopping integral t = leV 
(corresponding to a band- width of W = 8eV) , electrode- 
molecule coupling is g — O.lcV, Temperature is room 
temperature, T = 300K and the electrodes are kept at 
half filling. 

The first thing to be noted is that whenever a molec- 
ular orbital crosses the Fermi energy, the molecular oc- 
cupation changes by one, and the thermopower changes 
sign. This is in accord with the known results from 
the Landauer formula for thermopower, and is due to 
the change from electron-dominated to hole-dominated 
transport every time the Fermi level is crossed (which 



also corresponds to a transmission resonance). The sec- 
ond important feature is that, in contrast to t he result 
expected from the regular Landauer formulaP^, the See- 
beck coefficient exhibits strong variations with gate volt- 
age. 

Before we proceed to discussing the origin of these vari- 
ations, it is useful to compare the properties of the See- 
beck coefficient as obtained from the OQS method to 
those known from the Landauer formalism. In Fig. [IJb) 
the Seebeck coefficient as a function of gate voltage 
(same as in Fig. [TJa)) is plotted for different values of 
the coupling between the electrode and the molecule, 
g = 200, 100 and lOmeV. As expected from the Lan- 
dauer formalism 24 , there is little effect to the coupling 
on the magnitude of the Seebeck coefficient. In Fig. [TJc) 
the Seebeck coefficient is plotted as a function of tem- 
perature for various molecular chain lengths. An inho- 
mogeneous temperature dependence is found, originating 
from a crossover from coherent to incoherent transport, 
again in agreement with results obtained from the Lan- 
dauer formalism (e.gP). These results demonstrate that 
the main physical features of the thermopower which are 
present in the Landauer formalism, also appear within 
the framework of the OQS theory. 

We now turn to calculating the distribution of the 
thermo-voltage AV across the molecular junction. To 
obtain the distribution, we calculate the temperature- 
difference induced voltage AV — S(E)AT, taking the 
gate voltage of the molecule (i.e. the HOMO-Fermi en- 
ergy offset) £ to be a random variable, normally dis- 
tributed around Eq — 0.3eV with a width = 0.2eV. 
While these values are rather tdifferent then experimen- 
tal values (probably Eq is bigger in experiments) we point 
that we are aiming at qualitative similarity to experi- 
ment, to point the origin of the variations, rather than 
to analyze realistic junctions. To mimic the experiment, 
we also apply a small variation (two degrees Kelvin) to 
the electrode temperatures (note that in the experiments, 
a finite thermo-voltage distibution was observed even at 
AT sa 0,see inset of Fig. [2](b) , indicating the existence of 
a small temperature difference, probably due to Johnson 
noised This does not affect our results or conclusions). 

In Fig. [2^ a) the distribution of thermo-voltage is plot- 
ted for temperature differences AT = 0, 5, 10, 20, 30K. 
The distributions qualitatively resemble the experimen- 
tal distributions, exhibiting a broad double-peak struc- 
ture. To understand the origin of the distribution shapes, 
in the right panel of Fig. [2ja) the Seebeck coefficient is 
again plotted (solid line). Due to the strong sensitivity 
of S on E, even a relatively small variation in E (of the 
order of Te ~ 0.2eV) can include several maxima of S, 
giving rise to the different peaks in the distributions. 

To understand the origin of the sensitivity of S, in the 
right panel of Fig.[2ja) we plot the local density of states 
at the point of contact between the electrode and the 
molecular chain (dashed line) . One can see that the vari- 
ations of S(E) and the LDOS vary on the same energy 
scale. Note that these variations are a surface effect, and 
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are not due to the electrode level spacing (the average 
level spacing is ~ O.OleV). We thus conclude that it 
is the variations in the LDOS (which are on an energy 
scale much smaller than the energy difference between 
the molecular levels) that give rise to the sensitivity of S 
(we note that the total DOS of the electrodes is a much 
smoother function with no observed variations on these 
energy scales). We point that the LDOS variations are 
not a finite-size effect, but rather are a surface effect. 
Calculating the LDOS to systems as large as 200 x 200 
(where the level spacing is ~ 0.0005eV), we found that 
local variations in the LDOS persist with roughly the 
same energy scale of ~ O.leV (although their magnitude 
somewhat changes). 

On the other hand, when considering the thermopower 
as obtained by the Landauer formula, the local variations 
in the LDOS are not taken into account, or rather they 
are smeared by the electron self energy which is reflected 
through the level broadening Y in Eq. Q . The resulting 
S(E) is a much smoother function. As a first demon- 
stration of this effect, here we use the formula for trans- 
port and thermopower through a resonant level, Eq. [T] 
The positions of the levels and the level broadenings are 
obtained by fitting the data of Fig. [TJa). In the right 
panel of Fig. [2jt>) S(E) is plotted using Eq. p} and the 
parameters of the junction (i.e. resonances and widths 
fitted to Fig.[jja)). For comparison, S(E) obtained from 
OQS theory is plotted as a bright line in the background. 
The resulting Ay distributions are plotted in Fig. [2jb) 
for different AT. The distributions obtained using Eq. [2] 
are quite different from the experimental results (inset of 
Fig. |2pb)), both in terms of shape and in the lack of the 
double-peak structure. 



IV. GREEN'S FUNCTION ANALYSIS OF THE 
TOY MODEL 

Let us now examine the same model within the Green's 
function formalism. In this formalism, and for our non- 
interacting toy model, the transmission function t(E) is 



g 

given b-yiasiasj 



t(E) =Tr (Y L G r (E)Y R G a (E)) 



(9) 



where G r ' a (E) = {E — % =F iS) are the retarded and ad- 
vanced Green's functions, and Yl.r represent the level 
broadening due to the electrodes, typically a few hun- 
dreds m eVs^Z122l. j n t ne b as i s of atomic orbitals (as the 
Hamiltonian is written) Y^ r are diagonal matrices, with 
TL j n(n,n) = r if n is in the left or right edges of the 
electrodes, and zero otherwise. We take Y — 0.3eV in 
the numerical example below. £ is the self-energy, which 
in the non-interacting case is only due to the electrodes, 
and hence £ = Yl + Tr. 

Once the transmission is calculated, the thermopower 
can be calculated directly using Eq. (JT|). In Fig. [3ja) 
the transmission function t(E) and the its logarithmic 
derivative (proportional to the thermopower S(E)) of the 
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Figure 2: Thermo-voltage AV histograms for different tem- 
perature differences AT = 0, 5, 10, 20, 30K, obtained by tak- 
ing AV = S(E)AT with E normally distributed around 
E — 0.3eV with width Ye = 0.2eV. (a) Histograms obtained 
from the thermopower S(E) obtained using OQS theory (right 
panel solid line). The LDOS at the point of contact between 
the molecular chain and the electrode is also plotted (right 
panel solid line), (b) Same as (a), but for the Seebeck co- 
efficient S(E) obtained using the Landauer formula, which 
does not exhibit strong focal fluctuations. Inset: Experimen- 
tal histograms of AV for a DMTBDT molecular junctions, 
taken with permission from RefP. 



toy model are plotted (as a function of the gate energy) . 
For comparison, the local density of states at the point 
of contact between the wire and the electrode is plot- 
ted in the upper inset (the dashed line in the inset is 
the DOS when broadened by Y = 0.3eV). As seen, the 
oscillations in the LDOS are completely smeared in the 
transmission (plotted as a function of the gate energy). 
For comparison, the local density of states is plotted in 
the upper inset (the dashed line in the inset is the DOS 
when broadened by Y = 0.3eV). As seen, the oscillations 
in the LDOS are completely smeared in the transmission 
function and hence in the thermopower. The lower inset 
shows the transmission function on a log scale, verifying 
that the variations are really smeared out. To demon- 
strate that this smearing has a large effect on the ther- 
mopower, we consider an artificial transmission function 
f(E) = t{E) (1 + 0.1 cos(E/E ) with E = 0.02eV. This 
transmission exhibits oscillations on an energy scale Eq . 
However, they are chosen in such a way that they are 
hardly visible if the transmission is observed in its full 
scale, as seen from Fig. |3jb), where f{E) is plotted. In 
the inset the transmission is plotted on a log scale, and 
only there the oscillations are visible. On the other hand, 
the thermopower exhibits very strong oscillations of con- 
siderable size even if the oscillations are hardly observed 
in the transmission functions. This is due to the extreme 
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Figure 3: (a) Transmission function t(E) and thermopower 
S(E) as a function of energy offset for the tight-binding 
toy model. The transmission and the thermopower are 
perfectly smooth and disregard any variations due to the 
DOS. Upper inset: the local density of states at the point 
of contact between the wire and the electrode (sloid line), 
and the same LDOS broadened by V (dashed line), (b) 
same as (a) for the artificial transmission function f (E) = 
r(E) (1 + 0.1 cos(E/Eo). The oscillations in the transmission 
function are hardly visible, yet they induce large variations in 
the thermopower S(E). Inset: f (E) on a log-scale, only there 
the oscillations are visible. 



sensitivity of S(E) to local variations of t(E), reflecting 
its derivative structure. 



V. SUMMARY AND DISCUSSION 

In summary, we have analyzed the variations in the 
thermoelectric response of a metal-single molecule-metal 



junctions, based on detailed examination of experimen- 
tal results. The experimental results, namely the width 
of the thermo-voltage variations and the shape of the 
thermo- volt age distributions, cannot be accounted for by 
only assuming variations in the misalignment between 
the molecular HOMO and the electrode Fermi energy. 
Using the theory of open quantum systems we qualita- 
tively reproduce the experimental results, and show that 
they may originate in a combination of the level mis- 
alignment variations and variations of the local density 
of states at the point of contact between the molecule and 
the electrode, specifically the STM tip in typical molec- 
ular junction experiments. 

To put it differently, we found that in order to ex- 
plain the thermopower variations in molecular junctions 
the electronic transmission function cannot have a simple 
Lorentzian form, but rather should have a more compli- 
cated form which includes variations on a scale smaller 
than the HOMO-LUMO gap. Such variations may be in- 
duced by the variations in the LDOS at the tip-molecule 
point of contact, originating from the STM tip structure, 
impurities or trapped states, etc. The use of the Green' 
function formula smears out all variations on energy scale 
smaller than the level broadening (typically a few hun- 
dred meVs). This smearing has a very small effect on the 
transmission and hence the conductance of the molecular 
junction, but have a strong effect on the thermopower. 
Thus, caution needs to be taken when using the NEGF 
method for calculation of themopower. We note a recent 
paper^ arguing that the NEGF method has limitations 
even when calculating the conductance of a molecular 
junction (although the argument is different then pre- 
sented in this paper). 

The effects of STM tip structure on local transport 
measurements have been extensively discussed in the 
STM literature^!] and were recognized to have an im- 
portant contribution to the overall charge transport, via 
the tip LDOS. Our results imply that a similar impor- 
tance of STM tip structure appears in thermo-electric 
measurements. Consequently, a direct extraction of rel- 
evant parameters such as the HOMO-Fermi level offset 
from thermo-electric measurements requires taking the 
electronic properties of the metal- molecule interface into 
account. 

The author wishes to thank J. Malen for valuable com- 
ments on the manuscript. 
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